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ABSTRACT 



The true submersible, capable of extended operating periods without 
surfacing, has long been desired. Present day submarines must surface 
in order to charge batteries. Celestial navigation, therefore, is available 
as a by-product of normal operating procedures. With the advent of the 
submarine capable of extended periods of operation submerged, it will be 
necessary to provide a method of navigation that will not require the sub- 
marine to surface, if full advantage is to be taken of her other character- 
istics. 

The Instrumentation Laboratory, Massachusetts Institute of Technol- 
ogy, proposes, as a solution to the navigation problem, the application of 
inertial techniques employing gyros to provide the inertial reference and 
stabilized platforms to provide position indication. The problem is thus 
one of measuring the change in a direction, that of the true vertical, with 
respect to inertial space. 

Inertial techniques are being applied to air navigation. For these ap- 
plications, it is assumed that the carrying vehicle departs from a stable 
base, where a well defined vertical and azimuth direction are available 
for "lining-up" the system prior to departure. A somewhat different prob- 
lem exists for a submarine underway at sea. Here, there are no well 
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established references available; hence, the boundary conditions are 
largely unknown. 

A study of possible systems for applying inertial methods to sub- 
marine navigation is currently underway at the Instrumentation Labora- 
tory, Massachusetts Institute of Technology. A report (15) has been 
issued on the initial phase of the study, in which only undamped systems 
employing ideal components were considered. A self-settling system, 
however, must employ damping, which, together with component per- 
formance, will affect the accuracy of the system; hence, the complete 
study will consider these effects on system performance. Furthermore, 
the resulting system must have a period approaching the 84. 4 -minute 
Schuler-tuning condition if it is to be effective. This fact implies that 
the transient effects due to disturbances imposed upon the system will 
have a duration measured in terms of hours. It would therefore be de- 
sirable to have some means available for predicting the errors resulting 
from such disturbances. 

This thesis suggests that the mean-square error in position can be 
predicted for any particular time during the duration of the transient by 
using the power spectral density of the disturbance and a function defined 
as the transient transfer function, and that this method may be applied in 
the design stage to arrive at component performance specifications. 

For illustration, the equations for one proposed inertial system, us- 
ing nonideal components which in themselves generate uncertainties or 
disturbances, have been derived and the application of the proposed meth- 
od for evaluating the resulting errors is shown. 

It is regretted that lack of time and non-availability of equipment pre- 
vented a more thorough study of the problem. However, from the work 
completed, it appears unsafe to design a system on the basis of steady - 
state servoanalysis techniques where the solution of the problem during 
the duration of the transient is also of primary importance. This conclu- 
sion is based upon the fact that the response of the third-order system 
studied to a forcing input e^^^^ made an initial excursion which was up 
to several times as large as the steady-state response. Furthermore, 
the ratio of the initial excursion to the steady-state response increased 
with an increasing frequency, whereas the time at which the maximum ex- 
cursion occurred did not appear to be any well behaved function of the fre- 
quency. 
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OBJECT 



To investigate the possible application of a self- 
contained inertial system to the navigation of sub- 
merged vehicles. 
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CHAPTER 1 

DISCUSSION OF INERTIAL NAVIGATION 
INTRODUCTION 

The problem of navigation of a submerged vehicle is of interest 
in view of proposed submarine operating characteristics that enable 
extended submerged periods. These characteristics will make it nec- 
essary to navigate from the submerged vehicle for considerable per- 
iods of time, with some degree of accuracy, and without reference to 
external information. Solution of the problem of long-range all-weath- 
er aircraft navigation has been attempted by inertial techniques. An 
attempt to apply similar techniques to the navigation of underwater 
vehicles is only now being investigated. 

Six systems for inertial navigation have been proposed and studied 
in the undamped form, using ideal components. This thesis pro- 
poses to investigate the effects on the performance of one of these sys- 
tems when nonideal components, which introduce uncertainties such 
as gyro drift, integrator drift, etc. , are included in the system, and 
further, to extend the study to the damped system. 

INERTIAL NAVIGATION 

Navigation is the science of guidance of a vehicle over the Earth, 
over or through its waters, and through its atmosphere. This guid- 
ance involves a two -fold function of indication and control. Indication 
is the comparison of the indicated position of the vehicle with some 
reference fixed in the Earth. This comparison may be used to obtain 
the distance to a destination and the direction the vehicle must move 
to reach that destination. It has thus effectively produced the control 
signal. 
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Piloting, which is the simplest form of guidance, is ignored in 
this discussion, because it is assumed that no optical line of sight is 
available. The substitution of electro-magnetic beams for optical 
lines of sight, and the substitution of time intervals for distance in- 
formation, leave the essentials of piloting in the guidance problem. 

This guidance system is subject to interference, either natural or 
man made. The desire to minimize interference has stimulated the 
search for a self-contained guidance system operating on such natural 
earth phenomena as are effectively interference proof. 

Inertial navigation can be regarded as a precision dead-reckoning 
method in which measurements are made with respect to inertial space 
rather than with respect to the water surrounding the vehicle. Conven- 
tional celestial navigation is a form of partial navigation in which stars 
represent the inertial space references; but the requirement for sub- 
merged operation precludes the consideration of celestial systems. Full 
inertial systems use gyroscopes for inertial reference. 

Physical Factors Available for Guidance 

The natural phenomena associated with the Earth that offer possibi- 
lities for use with a self-contained guidance system include data from the 
medium surrounding its vehicle, the Earth’s electric field, the Earth’s 
magnetic field, and the Earth’s gravity field. References (13) and (1) 
discuss each of these possibilities in general and conclude that the Earth’s 
gravity field appears to be the best -suited source of data for a self- 
contained guidance system. 

The mass of the Earth has associated with it a gravitational field, 
which field is essentially radial to the center of the Earth and has no 
poles. The direction of this gravitational field could be readily deter- 
mined by a plumb bob if the Earth were non-rotating. Due to a basic 
law of physics, namely, the principle of equivalence in the general 
theory of relativity, that gravitational mass is equivalent to inertial mass, 
the daily rotation of the Earth produces a centrifugal acceleration that 
also affects the plumb bob. The vector resultant of the Earth’s gravita- 
tional field and the centrifugal force field is defined as the Earth’s gravity 
field, and the direction of this essentially radial gravity field is unique 
at any given point on the Earth (Fig. 1-1). 
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Fig. 1-1 Gravitational field intensity and gravity field Intensity (the vertical). 



The Earth itself is subject to the equivalence of gravitational and in- 
ertial effects. The associated centrifugal force field causes the Earth to 
bulge at the equator, giving it a nearly ellipsoidal or spheroidal shape. 

The figure that a fluid body with the mass distribution and daily rotation 
of the Earth would have is defined as the geoid. Attempts are made to rep- 
resent the geoid by an associated reference ellipsoid for the basis of 
geodetic measurements and map data. The surface of the geoid is rep- 
resented by mean sea level, which is everywhere an equipotential sur- 
face of the Earth's gravity field. Because the geoid does not have a 
consistently smooth surface, the vertical is not in general parallel to the 
normal to the reference ellipsoid at the same position. The angle between 
the vertical and the ellipsoid normal is called the deflection of the vertical . 
This deflection is generally less than 0.30 milliradian ; and generally, the 
largest marine deflections occur in regions of rapid changes in ocean depth. 

In general, variations in the deflection of the vertical can be consid- 
ered to be quasi-static at submarine speeds, and can be expected to have 
no dynamic effect on the operation of submarine inertial navigation systems. 

Basic Principles of Inertial Navigation 

The vertical, which is the direction of the specific force of gravity, 
is unique at any point on the Earth and is the basis of the astronomical 
position of that point. Astronomical position is a direction relative to the 
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Earth and has two components: 

1. Astronomical latitude, the angle between the equatorial plane and 
the vertical at the located point. This quantity is unique because the 
Earth’s polar axis has a unique spatial direction. 

2. Astronomical longitude, the angle between the plane through the 
Earth’s polar axis containing an arbitrary vertical (e. g, , that at Green- 
wich) and the plane through the polar axis containing the located vertical. 
This quantity itself is not unique because of its arbitrary reference, but 
difference of longitude is unique. 

Astronomical position data will not agree, in general, with geodetic 
position data (on which maps are based) due to the before mentioned de- 
flection of the vertical. However, position data determined by vertical 
indication and corresponding map data should agree in the case of small 
land masses whose positions have been plotted from astronomical data. 

The fundamental problem associated with inertial navigation is then 
the indication of the vertical inside a self-contained vehicle that is under- 
going additional linear accelerations. It is this problem which has been 
considered in detail in Reference (15). 

Theory of Vertical Indicators 

The simplest vertical indicator is the plumb bob or pendulum. A pen- 
dulous element is defined as a system free to rotate about a point, so that 
it assumes a preferred orientation under the influence of impressed ex- 
ternal forces. If the suspension -point is not accelerated with respect to 
the Earth, the pendulous element can provide a precise indication of the 
true vertical, i. e. , the direction of gravity. If the suspension-point is 
accelerated with respect to the Earth, the pendulous element tends to in- 
dicate the apparent vertical, i. e. , the direction of the resultant specific 
force. 

When the pivot of the pendulous element is accelerated, a reaction 
torque acts on the center of the mass of the element, producing a torque 
which causes the element to lag behind the pivot. To study the effects of 
an acceleration without introducing the effects of the mass of the element, 
it is convenient to discuss the motion in terms of force per unit mass, i. e. , 
in terms of the inertia reaction specific force associated with the accelera- 
tion. An inertia reaction specific force has the magnitude of its associated 
acceleration (though expressed in different units) and is opposite in direc- 
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tion. Figure 1-2 is a force diagram for an accelerated pendulous ele- 
ment in a gravitational field, indicating that the direction of the .result- . 
ant specific force, i. e. , the apparent vertical, is the vector resultant 
of the gravitational specific force and the inertia reaction specific 
force. 



neiu iniensicy) 

Pendulous Element Seeks to Align its Pivot — Center of Mass 
Separation with the Resultant Specific Force. 

Fig. 1-2 Specific forces acting on a pendulous element. 

To obtain the specific forces acting on a pendulous element during 
a navigational problem, it is necessary to know the associated accelera- 
tions, The linear acceleration of a point with respect to a reference 
coordinate system may comprise five components with respect to a 
coordinate system that is moving relative to the reference coordinate 
system. However, as far as marine navigation is concerned, only two 
components are significant. 

1. The tangential (horizontal) acceleration associated with the angu- 
lar acceleration of the true vertical as the vehicle accelerates over the 
Earth’s surface. 

2. The centripetal acceleration parallel to the projection of the 
Earth -radius on the equatorial plane, due to the Earth’s daily rotation; 
this, added to the gravitational field intensity, gives the gravity field in- 
tensity. 

When a vehicle moves over the nearly spherical earth, the vertical 
associated with its instantaneous position rotates approximately geocentri- 



Pivot — Center 




of Mass 
Separation 



Pivot 



Linear 
Acceleration 
of Center of 
Mass 



5 



cally with respect to the Earth. On a ship moving at constant speed along 
a great circle over the Earth, a pendulum would have to rotate with re- 
spect to the Earth at an effectively constant angular velocity in order to 
indicate the vertical continuously. Once the transient stages were over, 
no torque would be required to keep the pendulum along the vertical dur- 
ing this rotation. If the vehicle were to accelerate, by changing either 
speed or heading, a torque would be required to maintain the pendulum on 
the vertical, which in turn would be geocentrically accelerating. 

Schuler pointed out in 1923 that accurate indication of the vertical 
during periods of acceleration could be realized by suitable tuning of the 
natural period of the pendulum. Whenever the pivot of a pendulum is ac- 
celerated, the center of mass of the pendulum tends to lag behind the piv- 
ot with respect to inertial space. At the same time, this acceleration 
causes the true vertical, associated with the pivot, to accelerate geocen- 
trically with respect to inertial space. These considerations lead to the 
following observation: 

If a pendulum initially hangs vertically, it will remain along the ver- 
tical if its angular acceleration about its pivot equals the geocentric ang- 
ular acceleration of the vertical. (The accelerations are with respect to 
inertial space.) 

These two accelerations become equal for a distributed-mass pendu- 
lum when the ratio of the pivot -center of mass separation to the square 
of the radius of gyration of the pendulum equals the reciprocal of the radius 
of the Earth. When operation is in the Earth’s gravity field, and the pen- 
dulum is undamped, this condition gives a pendulum with a natural period 
of approximately 84 minutes. For a shorter -period pendulum the rotation- 
producing torque on the pendulous element is stronger, so the pendulum 
lags the true vertical. For a longer-period pendulum the torque is weaker 
and the pendulum leads the true vertical. 

An 84-minute period is very long compared to the periods of physical 
pendulums ordinarily encountered. In fact, it is improbable that a simple 
(concentrated-mass) pendulum or a distributed -mass pendulum could be 
constructed with this period, although the 84-minute gyropendulum just 
fails of being realizable. However, in the practical case, a condition anal- 
ogous to the Schuler -tuning of a pendulum can be realized in a particular 
kind of closed-loop system in which a short-period pendulum tracks the 
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the specific forces. A platform, controlled on the basis of the pendulum 
output data, then becomes an equivalent Schuler pendulum. 

From the foregoing, it is deduced that a vertical indicating system is 
conceivable which, while not a pendulum, utilizes the Schuler -tuning idea, 
so that a torque caused by the vehicle’s acceleration acting on the system 
will be absorbed by the inertia reaction torque required to keep a refer- 
ence line in the system on the vertical. 

Reference (15), from which the foregoing section was obtained, has 
elaborated upon this theory and set forth six representative inertial sys- 
tems for position indication. 

Fundamentals of Position Indication 

Astronomical position measurement involves the determination of the 
essentially geocentric angles between indicated and reference verticals. 

It is proposed to determine these angles by using a stabilized platform as 
an essential feature. The stabilization loop receives specific force data 
from accelerometers or pendulums, and delivers these data to gyro units 
and servo drives, which in turn drive the platform toward a horizontal 
position. The specific force data are modified and filtered to control the 
system dynamics. Part of this modification consists of Schuler -tuning 
the stabilization loop. Two basic methods of utilizing the indicated ver- 
tical are then available to obtain the vehicle’s essentially geocentric angle 
of travel. One, called here analytical integration, involves the integra- 
tion of accelerometer outputs or of gyro inputs, with the accelerometers 
and gyro units mounted on the stable platform. The other method, called 
here geometric integration, involves the direct comparison of the indicated 
vertical with a reference vertical. Both methods have been outlined^^^^ 
in terms of six representative systems utilizing ideal components and un- 
damped systems, namely: 

(1) Open -Chain Integration of the Angular Velocity of the Indicated 

Vertical. 

(2) Open-Chain Double Integration of the Angular Acceleration of 
the Indicated Vertical. 

(3) Direct Double Integration of Acceleration. 

(4) Electromagnetic Composition of Tangential and Vertical Com- 

ponents of Celestial Longitude Rate. 
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(5) Geometric Integration Using Angular Velocity of the Indicated 

Vertical with Respect to Inertial Space. 

(6) Geometric Integration Using Pre -aligned Inertial Reference. 

The first system will now considered in some detail using nonideal 

components in both the damped and undamped systems, with special em- 
phasis on the problem of determining a possible method for setting up 
performance specifications for components. It should be noted here that 
the normal methods of frequency analysis do not apply directly where the 
system does not reach steady state during its normal operating cycle, or 
where the transient period of the system is of primary interest. Such is 
the case in the operation of an inertial position indication system. 

Open-Chain Integration of the Angular Velocity of the Indicated Vertical 

Figure 1-3 is a functional diagram of a single -degree -of -freedom 
system for indicating the vertical. Figure 1-4 is a functional diagram of 
the essential operating components of a latitude indicator which uses an 
open-chain single integration of the angular velocity of the indicated ver- 
tical to arrive at the indicated latitude. Figure 1-4 corresponds to Fig. 
B-1 in Appendix B. Note that Fig. B-1 is primarily the basic vertical in- 
dicator of Figs. A -3 and 1-3, with the angular velocity integration tapped 
off the vertical indicator loop at the angular velocity signal that calls for 
processing the gyro unit, so that the position integrator output is the in- 
dicated change in latitude. The vertical indicator is undisturbed by the 
presence of this integrator, as far as the loop d 3 mamics are concerned. 
Hence, the same Schuler -tuning conditions applies to this system as ap- 
plied to the basic vertical indicator of Fig. A -3. 

A two -degree -of -freedom system constructed along these lines would 
require an x-axis accelerometer in addition to the y-axis accelerometer 
shown in Fig. 1-4, as well as x-axis components corresponding to the 
orientational control -signal generating system, gyro unit, and controlled 
member drive. In addition, the x-axis system requires that its angular 
velocity signal be modified by secant of latitude, in order to give longitude 
rate data to its position integrator. The position integrator in the x-axis 
system would indicate longitude difference with respect to some arbitrary 
reference. The maintenance of the y-and x-axes along a parallel of latitude 
and a meridian respectively must be handled by an auxiliary orientation 
system, e. g. , a gyrocompass. (1^) 
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1.3 Functional diagram showing interrelations among essential operating components of a 
single-degree -of freedom system for indicating the vertical when the gravity field vector lies in 

the plane normal to the Y-axis. 
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Fig. 1-4 Functional diagram showing interrelations among essential operating components of a 
latitude indicator, using an open-chain integration of the angular velocity of the indicated vertical. 



The gyro units ideally maintain their spin axes fixed with respect to 
inertial space. To maintain the controlled member on which they are 
mounted tangential to the Earth in its daily rotation therefore, the x-axis 
gyro must be precessed by a correction torque that keeps the platform up 
with Earth rate. The calculations of Appendix B consider this torque to 
be generated by a special precession current applied to the torque gener- 
ator of the gyro unit. 

In Appendix B, it is shown that, if the controlled member has the 
equation of motion of a Schuler -tuned pendulum, the correction to the in- 
dicated position (that is, the instantaneous difference between the true and 
indicated positions) has a similar equation of motion, except for the ap- 
pearance of certain initial conditions as additive constants of integration 
in the indicated position. Specifically, the correction to the indicated 
position at any time is the correction to the indicated position at the 
start of the problem, plus the present correction to the indicated vertical. 
The indicated position is thus Schuler-periodic with respect to boundary- 
matching terms associated with the vertical indication. The above descrip- 
tion holds directly with ideal components throughout the system. Position 
indication will be modified somewhat however, where uncertainties are 
present in the system, as is shown in Appendices A and B. 
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CHAPTER 2 



PROPOSED METHOD FOR ARRIVING AT COMPONENT SPECIFICATIONS 
UNCERTAINTIES IN COMPONENTS 




Fig. 2-1 Graphical definitions of inaccuracy, roughness, and error^for a 
quantity varying with time. 



Figure 2-1 is a graphical definition, by Draper of inaccuracy, 
roughness, and error for a quantity varying with time. The relationships 
shown are summarized below for convenience. 
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In the remainder of this thesis, when the error and correction in q 
are referred to, it will be taken to mean (E)Q(gm) and 
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The ideal components used in the inertial systems described in Re- 
port R-9 have been replaced, for the purpose of this study, with com- 
ponents which are no longer ideal in that the individual components now' 
are considered as introducing spurious signals into the system by their 
operation. 

For the purposes of this study, uncertainty in the output of an 
integrator has been treated as an equivalent spurious input to the inte- 
grator The specifications for an integrator in which the process 

has no specified maximum value may be written: 

(U)q(out)(lnt) ^ (2-6) 

That is, the specification is such as to require that the uncertainty of 
the output quantity of the integration process be less than a percentage 
of the magnitude of the integral of some specified value of the input 
quantity to the integrating system, over the interval of integration. In 
this connection, it must be assumed that the uncertainty of the integration 
process will not be truly random over the period of integration because 
of the fact that the integrator may be subject to inaccuracies, such as 
drift and frictional losses, w'hich do not lend themselves to correction; 
therefore, the average uncertainty of the integrating process will be 
different from zero over the period of integration. Similarly, the un- 
certainty in the output of an amplifier may not be a truly random pro- 
cess over a given interval because of amplifier drift, and the average 
uncertainty will therefore be different from zero over the interval. 

From the foregoing, it appears that specifications similar to those 
placed upon the integrator may be placed upon the amplifier, i.e., 
that the uncertainty in the output shall be less than some percentage of 
the magnitude of the output due to a specified value of the input quantity 
taken over a specified period of time. 



The foregoing may be represented mathematically as in Equation 
(2-7 ) where 
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It should be noted that, in both the integrator and the amplifier, the 
specification implies that the input signal has undergone an ideal inte- 
gration or amplification, and that superimposed upon the output is an in- 
accuracy which is due wholly to an uncertainty; specifically, that the pro- 
cesses of integration and amplification are carried through without error. 
Uncertainties in gyro units and controlled member drives, pendulum units, 
and accelerometer units as developed by the Instrumentation Laboratory 
are generally susceptible to separation into two categories: those due to 
uncertainties in electronic signals applied as inputs to the torque summing 
member or derived as outputs from the position of the torque summing 
member; and those due to torque uncertainties acting on the torque sum- 
ming member, such as mechanical unbalances, bearing friction, and 
undesirable reaction torques applied to the torque summing member by 
the microsyn signal units, etc. The subject has been treated extensively 
in several references, (10), (7), (6), and (16), and is of such scope that 
no attempt will be made to analyze it here, beyond the general observa- 
tion that the results of open-chain tests would seem to indicate that the 
powerful methods of statistical analysis might be necessary in order to 
determine the suitability of units of various designs for inertial system 
applications. This method, which is essentially that recommended in later 
sections, again would treat the units as ideal, with uncertainties superim- 
posed on the output of the ideal unit; hence, preliminary specifications may 
be laid down following principles to be developed later. 

In the previous section, uncertainties are defined, and the manner in 
which uncertainties may arise in various components are discussed. Now, 
it is necessary to consider the effect of uncertainties when the components 
from which they arise are employed in closed -loop systems. In particular, 
performance of components in a closed -loop system, as for instance the 
vertical indicating loop of Fig. 1-3, affects the dynamics of the loop, there- 
by determining the precision of indicating the vertical. Components which 
operate in open -chain fashion outside the loop, as for instance the position 
indicating integrator of Figs.B-1 and B-2, have direct additional effects on 
the precision of the system. Hence, the open-chain specification method , 
as proposed in previous paragraphs, as well as the method for arriving at 
performance specifications for closed-loop components, is of direct ap- 
plication. 
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By the assumption that is operated on by a particular component 
ideally, but that the component itself then contributes an uncertainty in the 
output which is superimposed upon the true output, we have in effect sepa- 
rated the total output ihto the sum of two outputs, one of which is the true 
output, and the other, the portion of the total output due to the uncertainty. 

If, then, it were possible to predict by some simple means the magnitude 
of the uncertainty and the way in which it varies with respect to some vari- 
able, say time, then it should also be possible to determine what effect the 
uncertainty would have upon the performance of the system. 

Unfortunately, the very term uncertainty implies the nature of the dif- 
ficulty. We cannot by any simple means determine the necessary informa- 
tion on the behavior of the uncertainties. However, by statistical study of 
the behavior of components under various operating conditions, it may be 
possible to determine the probable form of the uncertainties. In this re- 
port, it is considered that the results of such a study would be presented in 
the form of power spectral densities, which would then be used to enter the 
response curves for a given system, in order to determine the effect on sys- 
tem performance of a component with an output uncertainty of a certain spec- 
tral density. 

Since, in a position indication system, we are dealing with a dynamical 
system of second -order or higher, and, since our systems are lightly damped, 
the system may have normal modes or resonant frequencies that can be ex- 
cited by forcing functions of the same frequency. In particular, Wrigley has 
shown that an undamped vertical indication system of the type described 
later in this thesis is insensitive to a forcing functionwith an 84.4-minute period. 

The mathematical expression of the weighting function describes these 

k o t 

modes, which are found to be the terms of the form t eM into which the 
weighting function can be broken. In a complicated system, it may be very 
difficult to obtain an explicit mathematical statement for W(t), and, in de- 
sign problems, the amount of work involved in optimizing system parameters 
for the presence within the system of certain exciting frequencies may be ex- 
cessive. A further discussion of the weighting function is included in later 
sections. On the other hand, as shown in Appendix D, it may be a relatively 
simple matter, using machine methods, to determine the steady -state trans- 
fer function, Y(jw), the weighting function and the transient transfer func- 
tion, F(w,t), which is defined and discussed later; and from these to con- 
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struct charts from which system specifications can be prepared. 

The purpose of the following sections is to present in detail the 
manner in which such information may be utilized either in preparing 
specifications or in predicting the performance of the system. 

Again, it should be emphasized that the necessity for evaluating 
position information at all times the system is in operation requires 
a detailed consideration of the transient which persists for a consid- 
erable length of time because of the long natural period of the system, 
and the long system characteristic time. 



THE TRANSIENT TRANSFER FUNCTION 

Consider some arbitrary system performance function that may 
be represented 



q(out) ^ (p + V 

q(in) + a^p + 



(2-8) 



The transfer function Y(p) is frequently defined as the expression on 
the right with initial conditions zero, which is also the Laplace trans- 
form of the weighting function, where the weighting function is defined 
as the response of the system to a unit impulse with the system initially 
at rest. 

In the usual frequency response studies of a dynamical system, the 
steady-state response to a periodic function is the item of interest. In 
an inertial navigation system, however, because of the long system 
periods involved, the system may and probably will be placed in use be- 
fore the steady-state has been achieved. Furthermore, the effects of sys- 
tem disturbances are present as transients for such a length of time that 
a study of the transient response of the system to periodic forcing functions 
appears necessary. 

In order to study this transient response, we will define a new para- 
meter such that, if a system represented by a performance function of 
the form of Equation (2-8) is subjected to a forcing function of the form 
cos wt + j sin wt = e^^^, the response may be represented as follows: 



eiwt 



system 



F(w, t)e^^ 



Fig. 2-2. System performance function. 
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where, for the steady-state, Lim F(w, t) -►Y(w), the transfer function 

t —►CO 

for the system in which Y(w) will now be defined as the steady -state 
response of the system to the forcing function The F(w, t) then 

represents the transfer function of the system during the transient por- 
tion of the response, and is the parameter in which we are interested; 
moreover, it may be obtained by subjecting the system, initially at 
rest, to forcing inputs of the form e^^ for a series of frequencies and 
obtaining the time response of the system. | F(w, t)p thus obtained may 
be plotted vs. time for various frequencies. From these plots, cross 
plots of |f(w, t)| ^ for various times vs. frequency can also be prepared 
(Fig. 3-17). The manner of obtaining |f(w, t) | ^ is of interest and is 
covered in Appendix (D). 

Now, postulating that any disturbance to which the system is likely 
to be subjected can be assigned a power spectral density G(w),we may 
write for the steady -state: 

G(out)(w) =|Y(w)| 2 G(j„)(w). (2-9) 

Analogously, for the transient period, 

=|F(w,t)|2G(j„)(w). (2-10) 

Further, it can be shown that the ensemble mean square value of the out- 
put quantity at any time tj, for any arbitrary input with power spectral 
density given by G^j^j(w), is the integral of the product of the spectral 
density of the input and the square of the transient transfer function at 
time tn between the limits w = 0 to w = oo, that is, 

MS output at time(tj)^ |f(w, tj)| dw. (2-11) 

It should be pointed out here that the ensemble mean square may be 
defined as the mean of the squared end-point values of q(t) at time t^ for 
a number of trials; whereas the time mean square commonly used in com- 
munications applications is the time averaged mean square of q(t) over a 
specified time interval. Since we are interested in the output of our sys- 
tem at any specified time rather than the average output over some arbi- 
trary time interval, it follows that the ensemble mean square of the out- 
put is the quantity of interest. 

Now supposing |f(w, t)| ^ to be determined by some such means as that 
described in Appendix (D); a plot might then be constructed presenting 
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|f(w, t) I ^ vs. frequency, as shown in Fig. 3-18. It will be noted that 
by using this plot, and assuming that may be represented as a 

constant, the integration of Fig. 3-17 could be performed as the product 
of and the area under the curve of | F(w, t)| Or more general- 

ly, given plots of | F(w, t) | ^ and spectral density of a forcing function 
G(in)(w) vs. wfor any system under consideration, neither curve being 
necessarily representable analytically, where | F(w, t) | ^ is plotted for 
time tj and spectral density is time invar lent, the integration of Equation 
(2-11) may be performed by any of several standard methods of graphical 
integration. 

An examination of Fig. 3-18 reveals that the open-chain uncertainty 
specifications previously discussed are insufficient to specify component 
performance in a closed-loop system, since no account is taken of pos- 
sible frequency components present in the uncertainties. Fig. 3-18, for 
example, shows that frequency components in the vicinity of FR = 1 will 
have far greater effect than,for instance, the same magnitude of uncer- 
tainty with frequency components in the vicinity of FR = 4. However, by 
making certain assumptions, it is possible to use the information obtained 
from a study of |f(w, t) | ^ to arrive at preiiininary component specifications. 

If the uncertainty output of a component is assumed to be a purely ran- 
dom process, i. e. , the uncertainty has the characteristic of ’’white noise”, 
the power spectral density of the uncertainty is then constant over an in- 
finite band of frequencies from zero to infinity. However, | F(w, t) | ^ may 
approach zero if the system with which it is associated has an upper cut- 
off frequency, that is, input frequencies beyond a certain upper value of 
frequency will have no discernible effect on the system output. All the 
systems considered in this report have this upper frequency cut-off char- 
acteristic, i. e. , the systems behave as low-pass filters, and the assump- 
tion of ’’white noise” may therefore be of use. 

It should be noted, however, that an assumption of ’’white noise” of a 
given amplitude, is one of the worst which can be imposed upon a system, 
since the integration then involves all frequencies with equal weight. It 
seems reasonable to assume that the actual operating component will give 
rise to uncertainties whose frequencies will be predominantly centered 
within certain frequency bands, with negligible components outside these 
bands. If then the maximum spectral density in the predominant bands 
were to be applied to the system as ’’white noise” of the same amplitude, 
the resulting system performance would be worse than the actual by the 
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value of the integral outside the predominant frequency bands. 

The above reasoning suggests that preliminary specifications might 
be written on the basis of the open -chain performance specifications if 
it is assumed that the uncertainties thus obtained are the worst possible 
’’white noise” cases. Then, any improvement which might occur be- 
cause of more favorable frequency distributions may be considered as 
a welcome gratuitous improvement in system performance. 

System specifications for position indica.ting systems are generally 
written as requiring that indicated position be within specified limits of 
true position during specified time limits. The procedures of the fore- 
going yield information on the mean -square errors in system indication 
due to uncertainties in individual components. It can be shown that the 
sum of the magnitudes of these individual mean -square errors is equal 
to the mean-square error in the system output when due regard is given 
to the possibility that the individual terms may be complex. Hence, the 
procedure followed in the remainder of this report is to analyze the sys- 
tems on the basis that each component in turn is nonideal, the remaining 
components in the system being considered perfect. The information 
thus obtained should be useful in evaluating which components contribute 
the most to system error, and in determining lower marginal require- 
ments which must be met by components to be considered for use in the 
system. 

THE WEIGHTING FUNCTION 

The weighting function (W(t)) may be defined as the inverse trans- 
form of the term on the right of Equation (2-8). For further discussion 

of this subject, the reader is referred to any good text on servomechan- 
(91 (21 

isms. ' ' ' ' The weighting function for any system can be determined 
by subjecting the system to a unit impulse, i.e., 

W(t) = (PF)3y3t,^ (Ui)t, 

where 

(uj)t = unit impulse. 

Conversely, since the unit impulse is the derivative of the unit step, 
the weighting function is also the derivative of the response of the system 
to a unit step input. 

The weighting function has the characteristic that it relates the pre- 
vious inputs to the system to the present output by the integral 
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(2-12) 



^(out)W=7 V) (t-r) w(r) dr, 
0 - 



where 

r = t - tj . 

Hence, a knowledge of the weighting function enables one to form an esti- 
mate of the transient errors present in the output of the system following 
a forcing function which is applied over a given period of time and then 
goes to zero, such as the transient errors in an inertial system caused 
by maneuvers of the ship. This conclusion follows, if it is assumed that 
any given maneuver, occurring over a period of time, produces a forced 
error in the output of the system. If the ship then returns to a steady 
course and speed at time (t^), the forcing function is no longer being ap- 
plied to the system, but the error caused by the maneuver remains in 
the system at time (t) as a transient related to the forced error at time 
(tj) by the value of the integral in Equation (2-12). The above example 
is a specialized application of the general theory by which a series of 
impulse or step functions may be used to approximate any arbitrary in- 
put to the system, and thus the normal response of the system to the in- 
put may be determined from the relationship of Equation (2-12). 



20 



CHAPTER 3 



ANALYSIS OF VERTICAL INDICATORS 



GENERAL 

The differential performance equation for the position indicating sys- 
tem being considered in this thesis is derived in Appendix B. The differ- 
ential performance equations for the associated vertical indicators, both 
damped and undamped, are derived in Appendix A. Equation (B-7) shows 
that the correction to the indicated vertical, [(C) V]^|. enters directly as 

one of the terms in the correction to the indicated position, [(C)P]^^ 
Consequently, in the future, when references are made directly to the 
performance of the vertical indicator, the reader may infer that refer- 
ence to position indication is implied as well. 

In this chapter, we will consider first the undamped vertical indicator, 
and finally the damped vertical indicator. The concept of a transient trans- 
fer function as introduced in Chapter 2 will be used as an aid in evaluating 
the effects of certain uncertainties upon the vertical indication. Since the 
results of any uncertainty spectral -density analyses on the various com- 
ponents have not been available to the authors, if indeed such studies ex- 
ist in usable form, only qualitative estimates will be made as to the rela- 
tive performances required of the various components. 

UNDAMPED VERTICAL INDICATOR 

Figure 3-1 presents in summary form the equations of |F(w,t)| ^ for the 
undamped vertical indicator. Figure 3-2 is a semi-log plot of Equations 
(3-5) and (3-6) vs. frequency ratio for a time of one-half hour. Photo- 
graphs of transient transfer functions of damped vertical indicators are 
shown in Fig. 3-3. These damped systems are to be discussed later in 
this chapter. It will be recalled that the mean -square errors obtained in 
the proposed method are the ensemble mean squares, and consequently rep- 
resent the errors at a particular instant of time. Further, by the previous 
definitions, the time interval is measured from the instant the forcing 
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OlTeo the perforuAce equation of a lyatem for which the retponee to a fopcln* Input of the for* e^*f* !• deelred, and recalling that e^*f* 
m eoe W^t * J eln W^t, the reeponie of the lyeten to e^^f^ le glTen by the conplex eun of the eei>*rate reaponiee of the tjretea to eoe W^t and iln H^t. 
We can applythe foregoing to the undamped vertical Indicator with the differential equation (Bq. A^32), which nay be written 



<l(ln)0 * ^ '^(in)l. (*- 32 ) 

where eo* ^(ln)2 ** appropriate quantltlee appearing In the tero 

and flret order derlvatlvei reepectlvely of equation (A>-32) In Appendix A. 

Denote by )(^)> reeponee of the eyetea to where the eubecrlpt (n) la the order of the dorlvatlve aiaoelated with 

^(ln)(n)' ^(n)(aln)^^^' reeponee of the eyeten to “ eln W^t. Then talcing Initial condltlone equal to eero, the eolation for Xq- 

(Ap>32) to the forcing function e^^f^ may be written ae the eum of the four teme below. 



coa W^t - coa 






'(o)(co.)(*> “ ■ 


*nX " V 


nl 


(3-1) 


J'(o)(aln)<*J “ 


^ ( *nl - V 


r) (ala Wjt - am W^t) 


(3-2) 


^l)(co.)(»> " 


(iTnij) 


(W, am Wjt - W^ am w^t) 


(3-3) 


^'(l)(ala)^*^ = 




( - coa Wjt ♦ coa W^t) 


(3-4) 



Beealllng the definition of jT(w,t)| • Chapter 2, and malclng appropriate eubetltutlona following the method of Appendix D, We obtain 

(3-5) 

l»(w.t)l*ij - IJ '(l)(.ln)‘l* - f V - **1 ’'nl‘5“ * <«* W " ^-6) 

X HM 

In an analagoue nanner» the tine eolatlone for the aero order tern of the damped vertical Indicator of Xquatlon (A*38) nay be written 
1 ^ _t_ I ^ 

'(,)(«..)<*>- )-j . a )■ , (1 . • *{[(^r * V] K: - v>" ’ * »")• v]} 

e , ^ iyi T5 =r^* a-^^^^ne* eln (W Jl - (D»)* ♦ *3)* 

'il- (DB)a [[(5^ - (DB)W^ J ♦ W^» (1 - (DB)*)] [[W^J (2(DB)»- 1) ^ 4 (DH)» (1 - (DH)a)]J 

' [v ♦(c^)l(o”>».. - cij)' »•.: ('- m>*5r *!%)■ ♦ v] [i».; - v>’ • * <“>• v]} 



V t j 1 1/2 -(ra)Wj^,t 

"na^l- (“)* |[(^- (DB)W^J“ e Wj (1 - (DB)a)][cWj^; (2(DH)»-1) e W,*]* e 4 (DR)» (1 - (DH)a)]J * •*» ( ne " < ^ 2"*3^ * 



-1 -1 2(I«)W W, 

where A^ = - tan * CTjj - tan * - y a . y a 






na ■ "f 

W^aVl - (DB)a j 2W^^; (DB)Vl - (DB)« 

XT^ - W„5 (2(DH)*.l) Vwp ■ 






The correepondlng |T(w,t] for Squatlona (3-7) and (3-8) la not Incliided here, but the reader can readily verify that the reeultlng function 
le aonawhat conplex, and that |r(w,t)|^^^ la net a particularly well behaved function of the forcing frequency. 



Fig. 3-1 Summary of time solution for vertical indicators. 
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lF(w,t)l 




Fig. 3-2 Transient transfer function squared of undamped vertical indicator 
for time t = 1/2 hour. 
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Fig. 3-3 Typical oscilloscope presentation of transient transfer function of damped vertical 

indicator. 








function is applied to the instant in time under consideration. 

Also, although strictly speaking the frequency form of the transfer func- 
tion is not defined for an unstable system the frequency form of the 
transient transfer function does exist by the definition in Chapter 2. 

The area under | F(w,t)| plotted on linear coordinates would rep- 
resent the integral jT |F(w,t)| dw with = 0, 0<w< oo, and 

is less than the area under the curve for |F(w,t)|^j plotted in the same 
fashion. Hence, for "white noise" for which G^.j^j(w) is constant for all 
frequencies, the term containing the first derivative produces the greater 
error in vertical indication. 

The first derivative term is associated with the uncertainty in gyro 
drift and with the inaccuracy in the speed compensation current. Con- 
versely, the uncertainties in the accelerometer unit and the integrator 
within the vertical indication loop are associated with the zero -order 
term. Therefore, the performance of the gyro unit and the measurement 
of water speed will probably be more critical than the performance of the 
accelerometer unit and the integrator. This conclusion will be especially 
true if the spectral densities of the uncertainties have a wide band width 
such that the dips in | F(w,t) | ^ for the first-order term do not have an ap- 
preciable effect in reducing the associated errors. 

The one -half-hour time was chosen as representative of possible sit- 
uations in which a damped system would not be used. These situations 
would be those in which the period of operation would be short with re- 
spect to the time required to damp out transient disturbances, and hence 
no advantage would be gained from the greater complexity of the damped 
system. 

DAMPED VERTICAL INDICATOR 

Next, we consider the damped vertical and position indication system. 
The parameters used are those that provide the optimum position indica- 
tion with respect to the angular velocity of the true vertical. The result- 
ing parameters are summarized in Table C-I. 

In Reference (8), the optimization was done on the opposite basis from 
that used in Appendix C; specifically, the acceleration error was eliminated 
and the jerk-error minimized. It is interesting to note that the attenuation 
of errors due to accelerations under the conditions of Reference (8) pro- 
ceeds at 6 decibels per octave with increasing frequency, whereas 
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by using the present methods, acceleration errors are attenuated at 18 
decibels per octave. 

Figures 3-4 through 3-15 present, respectively, the steady-state 
sinusoidal responses of the vertical indicator, the step responses, and 
the impulse responses or weighting functions, for DR = 0. 1, 0.3, and 
0. 5. As defined in Chapter 2, | F(w, t)|is an analogous function to|Y(w)|, 
the transfer function; and hence the area under the curve of|Y(w)| ^ has 
the same utility for the steady-state as the area under the|F(w, t)| ^ has 
for the transient period. Hence, from a study of the steady-state re- 
sponse curves and a comparison with Equation (A -38), we may conclude, 
by reasoning similar to that applied in the case of the undamped system, 
that the uncertainties associated with the gyro unit and water-speed com- 
pensation will probably be most critical. Proceeding in a similar man- 
ner, we may arrive at the following listing, which is essentially a list- 
ing of components in direct order of the excellence of performance re- 
quired, but without regard to the sensitivities which may be involved. 

(1) (a) Gyro unit 

(b) Water -speed compensation 

(2) Accelerometer unit 

(3) (a) Direct-channel amplifier 

(b) Earth -rate compensation 

(c) Feedback amplifier 

(d) High-frequency integrator 

(4) Low-frequency integrator 

It should be noted that, in the damped vertical indication system, six 
sensitivities appear. These six sensitivities are related in Derivation Sum- 
mary A-H by only four equations. Thus, we may impose two more condi- 
tions on the sensitivities, such that we arrive at six independent equations 
in the six sensitivities. It is apparent, then,that we may be afforded some 
leeway in adjusting sensitivities such as minimizing the errors due to some 
particularly objectionable component in the system. 

No mention has been made in this chapter of the effect on position in- 
dication of uncertainties in the open-chain integrator. However, as shown 
in Appendix B, the mean-square error in position due to the open-chain in- 
tegration may be developed in a manner analogous to that applied above. 

The remainder of this chapter is devoted to a discussion of the methods 
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Fig. 3-4 Steady-state sinusoidal response of three subsystems of damped 
vertical indicator, DR = 0.1. 
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Fig, 3-5 Steady-state sinusoidal response of three subsystems of damped 
vertical indicator, DR = 0.3. 
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Fig. 3-6 Steady-state sinusoidal response of three subsystems of damped 
vertical indicator, DR = 0.5. 
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Fig. 3-7 Steady- state sinusoidal response of two subsystems of damped 
vertical indicator, DR = 0.1. 
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Fig. 3-8 Steady-state sinusoidal response of two subsystems of damped 
vertical indicator, DR = 0.3. 
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Fig. 3 -9 Steady state sinusoidal response of two subsystems of damped 
vertical Indicator, DR = 0.5. 
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Fig. 3-10 Step response of damped vertical indicator, DR = 0.1. 
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Fig. 3-11 Impulse response of damped vertical indicator (derivative of 
step response), DR = 0.1. 
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Fig. 3-12 Step response of damped vertical indicator, DR = 0.3. 
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Fig. 3-13 Impulse response of damped vertical indicator, (derivative of 
step response), DR = 0.3. 
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Fig. 3-14 Step response of damped vertical indicator, DR = 0.5. 
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Fig. 3-15 Impulse response of damped vertical indicator, (derivative of 
step response), DR = 0.5. 
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used and the results obtained in studying the transient transfer func- 
tion for the damped vertical indicator. 

For the study of the damped vertical indicator, it was attempted 
to apply the method described in Appendix D, using a high-speed, gen- 
eral purpose analog computer. 

Two identical systems (Fig. 3- 16a) were placed on the simulator 
and the forcing functions cos w^t and sin w^t were applied, one to each 
system. The output of one system was then taken to the x-axis of an 
oscilloscope. The resulting. Lissajous pattern represented a plot of 
the sine response of the system vs. the cosine response, from which 
it was proposed to obtain the transient transfer function, |f(w, t)| . 

It was intended that the sine and cosine functions be developed in 
an harmonic oscillator (Fig. 3-16b) as a part of the set-up. The oscil- 
lator would then oscillate at the frequency determined by the w^ setting 
when excited by a step input from the computer master generator. 

Theoretically, such an oscillator would be possible using an un- 
damped second-order system, and the sine and cosine outputs plotted 
against each other would produce a circle on the oscilloscope. However, 
due to nonlinearities or other imperfections in the computer, it was 
found that the amplitude of the oscillations of both terms tended to change 
in magnitude. Consequently, it was found necessary to introduce either 
positive or negative damping into the system, depending upon the fre- 
quency, to maintain a constant amplitude. 

The damping, in turn, affected the sine and cosine outputs differ- 
ently; hence, it was necessary to include the ’’scope-attenuator” in each 
system to adjust the resulting ellipse to a circle. The use of an addi- 
tional attenuator in this position will not introduce distortions in the sys- 
tem responses, since ail signals are attenuated the same amount for a 
given setting. 

The ’’input gain control” was put into the input to the harmonic os- 
cillator to provide a calibrated gain control in order to increase or de- 
crease the system responses to usable size. In tests, changes in set- 
ting of this control seemed to produce no distortion in the circular dis- 
play of the sine vs. cosine pattern, and hence it may be assumed that 
the amplitudes of both were changed equally. 
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to K, 



to K, 





(b) Harmonic oscillator. 



Fig. 3-16 Block diagrams of high-speed analog computer. 



40 






The remainder of the set-up is clearly apparent, the remaining gains 
and attenuators being applied merely to provide more reasonable potentio- 
meter settings. 

In calibrating the oscilloscopes, one of which was used as a monitor - 
scope, and the other as a photo-scope, the necessary damping ratio was 
first set into the harmonic oscillator. Then, with and Kg set to con- 
venient values, the oscilloscope amplifiers were used to set in a two-inch 
diameter circle on the monitor -scope. Then, with reduced, and all 
other settings maintained constant, a two-inch circle was adjusted on the 
photo-scope using the oscilloscope amplifiers. The shift in in the 
calibrations was made so as to take advantage of the linear amplifiers on 
the oscilloscope to amplify small pictures, rather than alternatively ap- 
plying large forcing signals with the attendant danger of distortion due to 
over -loading the computer. 

After calibration, the oscilloscope amplifiers were left undisturbed. 

For each frequency thereafter, a two-inch circular plot of the sine vs. the 
cosine was obtained on the monitor -scope by manipulation of the damping 
ratio in the oscillator and the scope attenuators. Kg. The input gain was 
set at a convenient value for the above operation at each frequency and 
recorded, as was the gain used for photographing each response, in order 
that the responses might be reduced to a common reference. Calibration 
was completed by recording the two-inch, or unit-reference circle photo- 
graphically, in order to provide a direct reference for reducing the data. 
Using the above methods, photographs of the Lissajous pattern responses 
of each of the sub-systems were made at various frequencies. It was ap- 
parent, however, that the highest forcing frequencies which might be ap- 
plied corresponded to a frequency ratio of about five, because of the large 
attenuation in the system at higher frequencies and the consequent limita- 
ations imposed by the noise level in the computer. Thus the data could not 
be extended to the significant frequency region for which FR > 16. The data 
were then reduced from the Lissajous patterns, of which those shown in Fig. 
3-3 are representative. The time dots on these patterns signify 4-hour 
intervals. An examination proved these data to be unreliable, since steady- 
state amplitudes thus obtained failed to agree with steady-state amplitudes 
computed by other means. As a consequence, no effort will be made to 



41 



present any quantitative results. 

Alternatively, a qualitative summary of the results of a study of the 
Lissajous patterns from sub-system to sub-system and damping ratio 
to damping ratio is presented below. The points in the summary occur 
conijistently and are considered significant in the transient study. 

(a) An excursion of the response occurs that is greater than the 
steady-state response. In addition, the initial excursion may or may not 
be oscillatory. 

(b) The response then returns toward the steady-state in an oscil- 
latory manner, the frequency of the oscillations apparently increasing 
with increasing frequency. 

(c) The ratio of the initial excursion to the steady -state response 
apparently increases with an increase in frequency. 

(d) The time in hours at which the maximum excursion occurs is 
apparently a function of both damping ratio and frequency for any parti- 
cular system, and, for the systems studied, varied from six to forty 
hours. 

(e) The time of reaching the maximum excursion does not appear to 
be an orderly function of frequency for a particular damping ratio and sys- 
tem. That is, for increasing frequencies, the time at which the maximum 
excursion occurs may first show an increase and then a decrease. From 
the data obtained, no reliable conclusion seems possible as to the effect 

of damping ratio on the time for reaching the maximum excursion. 

(f) The nature of the responses seems to indicate the inadvisability 
of designing a system purely on the basis of conventional steady-state 
analysis techniques, if the response during the transient period is of in- 
terest. 

The previous statements are not considered by the authors to be a 
condemnation of the method. It is suggested that additional work, which 
could not be undertaken by the authors because of a combination of lack 
of time and non-availability of equipment, might yield interesting results. 
Conservatively, these results would show low-frequency behaviors of the 
sub-systems for the vertical indicator. It would then remain to determine 
the relationship, if any, that exists between the behavior at high and low 
frequencies. 

A method proposed by Dr. J. H. Laning, Jr. , of the Instrumentation 
Laboratory, Massachusetts Institute of Technology, would approximate 
the high-frequency transient transfer function using the response of the 
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system to impulses. This method will be investigated at some la*^er date. 

Another method, suggested by Dr. Sidney Lees, also of the Instru- 
mentation Laboratory, would make use of a special applicatic of the roots - 
locus method of servoanalysis to reduce the amount of labor involved in 
calculating a point by point response of the transient. 

Still another method, for which equipment was not available, would 
again use the analog computer. Here, the transient transfer function 
would be presented as a function of time by squaring and adding the re- 
sponses of the two systems, taking the square root of the sum, and ap- 
plying the resulting signal to one axis of an oscilloscope, the other of 
which is the time axis. 

The remaining material in this chapter is presented not with the idea 
of showing any quantitative information, but by using the data as obtained 
to show the manner in which the theory might be applied. 

Figure 3-17 is a linear plot of | F(w, t) | ^ obtained by the methods of 
Appendix D. It is from this plot that it is proposed that mean -square er- 
rors be determined by the integral of Equation (2-11). 

MS output = y“|F(w,tj)| dw. (2-11) 

1 0 ' 

Figure 3-18 is the more familiar semi-log plot of Fig. 3-17, on which 
is also plotted, but not to scale, the steady-state transfer function for com- 
parison. Thus, it is seen that the apparent spikes in Fig. 3-17 are just 
such as would occur if the transfer function squared were itself plotted to 
a linear scale. The shape of the curve at four hours is due to the fact 
that the maximum initial excursion has not been reached, and further, 
the particular system under consideration is one for which the time to 
reach the initial peak decreases with increasing frequency over the range 
considered. It is interesting to note that the remaining curves for 12, 24, 
and 36 hours generally follow the shape of the transfer function. 

Figures 3-19 through 3-21 are the plots of the areas under the curves 
of Fig. 3-17 for the five sub-systems defined by the terms on the right of 
the simplified form of Equation (A -38), shown below. 
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Fig. 3-19 Ensemble mean square errors in damped vertical indicator, 

DR = 0.1. 
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Time (Hours) 



Fig. 3-20 Ensemble mean square errors In damped vertical indicator, 

DR = 0.3. 
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Fig, 3-21 Ensemble mean square errors in damped vertical indicator, 

DR = 0.5. 
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Cy) [p] q(in)|j'" ®(<l(in)> Cy) [p]^(in)j'' 



+ S, 



^^(in)’‘'V^ [P] *1(111)2 



+ 




+ S, 



( + a2P) 



(A -38) 



(*l(in)'‘^V^ 



[P] *l(in)(2, 1)’ 



where [P] = characteristic differential equation for 
the vertical indicator. 



Following the definition of Equation (2-11), the plots of the areas 
may be assigned the special significance of the mean -square error in 
vertical indication at each time, t = 4, 12, 24, and 36 hours, for an 
input with power spectral density such that = 1 for 1S(FR)<4. 75, 

and is everywhere else equal to zero. 

The reader will note the similarity between the shapes of the error 
curves for corresponding systems at different damping ratios, and fur- 
ther, the apparent trend in the relative location of the curves at each 
damping ratio. However, considering the accuracy of the data involved, 
it is doubtful that more than a tentative conclusion can be drawn as to the 
relative importance of the sub-systems as error sources until further 
work is done along this line. 
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APPENDIX A 



VERTICAL INDICATION 



BASIC PRINCIPLES 

The following principles of vertical indication were set forth in Ref- 
er ence( 15), and are repeated here for completeness. Figure A-1 is a 
sketch of a three -axis vertical indicator. 

The vector angle representing the correction to the indicated verti- 
cal is expressed mathematically as 



where the quantities are defined as in Fig. A -2. 

Indication with Respect to Earth Axes — Accelerometer Tracking 

The rate of change of this angle with respect to indicated position 
axes, i. e. , in a set of coordinates that consists of the indicated verti- 
cal and two other mutually orthogonal axes arbitrarily oriented in the 
indicated horizontal plane, is 




(A-l) 



[{C)V](t i)p_ - W(gv)t -W(Ev)i ~’^(EP)Vj ^ , (A -2) 



where 




rate of change of indicated vertical cor- 
rection with respect to indicated position 
axes. 




angular velocity of true vertical with re- 
spect to the Earth, 




angular velocity of indicated vertical with 
respect to the Earth, 



^(EP)Vj 



angular velocity of indicated position axes 
with respect to the Earth about the indicated 
vertical . 
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X-AXIS ACCELEROMETER 
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Fig. A-1 Line schematic diagram showing essential geometrical parts of a complete vertical 
indicating system using two single-axis accelerometers and three single -degree -of -freedom 



VECTOR FORM 

0P_[(C)V]„,„ 



RESULTANT 
SPECIFIC FORCE 



UNIT VECTOR ALONG 
INDICATED VERTICAL 



UNIT VECTOR ALONG 
TRUE VERTICAL 

T. 




^EV)t 
ANGULAR VELOCITY 
OF TRUE VERTICAL 
WITH RESPECT 
TO EARTH 



APPARENT 

VERTICAL 



INDICATED 

VERTICAL 



(t,i) 

CORRECTION TO INDICATED 
VERTICAL WITH RESPECT TO 
TRUE VERTICAL. (Angle through 
which Indicated Vertical must 
be rotated to align it with 
True Vertical.) 



AS A VECTOR ANGLE 



[(c)v](t,i) = X i 



''t sinC 



t,i 






SINCE [(C)v](^^.j IS SMALL 



Fig. A -2 Geometrical relationships associated with the vector angle form of the correction to 
the indicated vertical with respect to the true vertical. 
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The magnitude of the indicated horizontal component of the resultant 
specific forces, i.e. , the quantity sensed by an accelerometer, is ex- 
pressed mathematically in Equation (A -3). 

^®^\res)Hj = j (®^(res)Hj | ° | Vj ^ ^®^(res)[ ’ (A -3) 



where 






indicated horizontal component of resultant speci- 
fic force, 

denotes magnitude of the enclosed vector quantity. 



Let the system operate such that the angular velocity of the indicated 
vertical with respect to the Earth is a function of the magnitude of the in- 
dicated horizontal component of the resultant specific force and directed 
so as to rotate the indicated vertical toward the resultant specific force 
vector, namely 



''^(EV)i ■ W) I ^Vj ^ (®^)(res) ’ 



(A -4) 



where 

(PF)^yjjj£ = vertical -indicating performance function for a 
specific force input and an angular velocity out- 
put. 

Substitution of Equation (A -4) into Equation (A -2), plus the relation- 
ships for the resultant specific force from Derivation 1, gives 



[(C)V](t,i)p = ”(P^)(Vi){f,W) 



ly X 



I g * (Wgp) 



EP^E ^ ^E 



- Wep X (Wep X Re) - 2 Wje X (Wep X Rg)] - W(Ep)Vj X [(C)V](J^ y . 

(A -5) 

The specific force components in Equation (A -5) can be conveniently 
resolved into components parallel to the true vertical and perpendicular 
to the true vertical (the arbitrary nature of axes in the horizontal plane 
requires no further subdivision into components), namely. 



g - ^ » 



^E " “ ^E » 
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'''^(EP)Vj^ ^(EV)t • 



Substitution of these relationships into Equation (A -5) gives the specific 
force in terms of components along the true vertical, along the angular 
velocity vector of the true vertical (which must be perpendicular to the 
true vertical) and perpendicular to both of these, namely. 



Vj ~ ^E^ (EV)t 



[(C)'V^j^ 1)P. = ''^(EV)t “ W)^Vj ^ 

- 2Re(Wie^ • "^(EV)p] ^(EV)t^E[^(EP)Vj ^^^lE • 
+ [(Wep) e X V^^**e|"^(EP)Vj ^ • 



(A -6) 



Carrying out the vector operation indicated in Equation (A -6), re- 
grouping terms and omitting the components parallel to the indicated 
vertical as not being applicable to the indication of the vertical gives 



[(C)V](J J)p_ + (®^)(res)Vj " ^(EV)t 

-(PF)(Vi)(f,W)^E^'''^(EV)PPj “ (^^)(Vi)(f,W)^Ef''''(EP)Vj 

+ 2(Wjg.Iy^) ] (ly^ X W^gyjj) - W^gpjy_ X [(C)V](j^ Jj , 
where 

(®«(res) Vj = 6 ~ ^E^^(EV)t ~ ^{EV)t^ ’ -8) 

= component of resultant specific force along the true 
vertical . 



Let be the following operator: 

f ' 

where denotes integration of a function of the indicated axes with 

respect to the time. Substitution of Equation (A -9) into Equation (A -7) 
gives Equation (A -10): 
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[(C)V ](t^ i)p^ + ^ (®f){res) Vj 1) = ^(EV)t 

-^(vi)(f,^/p^ ^E<^(EV)t)pj**‘ ■Vi)(^,^^ /p. *’e[^(EP)Vj 

+ 2(Wjg . Ty^)](Iy^ X W(Ey)j) dt - W^gpjy. X .) . (A -10) 

Differentiating (A -10) with respect to the time gives the differential per- 
formance equation in terms of functions of the indicated axes: 

[(C)V](t J)p^ + (®*)(res)Vj 1) 

= (1 -S(yi)(f^-^ ^E^^^(EV)t)p. “ ®(vi)(f,W) ^E [^(EP)Vj 

+ 2(Wie. Iy^)](iy. X W(j,y,j) - (W(Ep)y. X i)p.) , (A-U) 

where it is assumed that variations in R^, are slow enough to be treated 
as quasi-static and variations in (sf)jj.ggjy are either slow enough to be 
quasi-static or fast enough to be ineffective. 

For Schuler -tuning let 

®(Vl)(f,W) = ’ 



where R^ is a set value made as nearly equal as possible to the magnitude 
of Rg. If it is then assumed that R^ = R^,, Equation (A -11) becomes 

»» (®^)(res)V _ 

[(C)V](t^ i)Pj + ^ 1) " “ [^(EP)Vj ^ ^''^(lE) • 

X W(Ey)t) - (W(Ep)y^ X [(C)V](j^ i)) p^ . (A-12) 

The terms W^gpj„ (i^ x W^gy^j) and 2 (Wjg. Ty )(Iy x W^gyj^) are 

(®^(res)Vj ' 

small in comparison with , and can be neglected. This can be 

■p 

shown as follows: s 



^(EP)Vj> 



the magnitude of the rate of change of the meridian about 



the true vertical with respect to the Earth, is given by 
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sin (Co) tan (.Lat) , 



(A -13) 



W 



(EP)V. 



where 



Vg = Velocity of base , 

Co = course of basej 
Lat = latitude , 



W, 



(EV)t 






(A -14) 



and 



W, 



(EP)Vj ^ "^(EV)t 



~''^(EP)Vj ''^(EV)t, 



(A -15) 



|2(Wje. ly^ (iy^ ^''^(EV)p|~ ^^(EV)t^IE 
Assume that the base is traveling due east at 10 knots at 70” latitude: 



'^(EP)Vj ''^(EV)t 
(®^(res)Vj 



^s''^{EP)Vj ^(EV)t 



(sf) 



= 00.02 min. of arc^ (A-16) 



(res)V. 



2Wie sin (Lat) W^gy^). ^ ^lE ^(EV)t ^ qq 35 min.of 

^®^\res)Vj ^®^^(res)Vj ' arc. 



(A -17) 



Neglecting the above terms, Equation (A-12) reduces to 



(sf) 



[{C)V](t t)p. + 



(res)V. 



[(C)V](tj) = -(W(Ep)y.X[(C)V](^i)). 



(A -18) 



The term (W^ep)V- ^ i)) represents an orientation transfer, not 

a magnitude effect, and represents a secondary coupling between the x and 
and y systems. Since this will only produce a ’’spiraling in” effect as the 
system approaches its null position, it will be ignored in subsequent dis- 
cussions. 
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Equation (A -18) now reduces to 

^®^^(res)V. 

— i ^ =0- 



(A -19) 



INSTRUMENTATION OF A SINGLE -AXIS OF THE VERTICAL INDICATOR 
OF FIGURE A-l(UNDAMPED) 

Definitions of Symbols Used in the Following Derivations 



‘"EB 

^IB 

Pi 

Pt 

[(C)P](t,i)o 

[(C)V](t^,) 

f(C)V](t,r)o 

®( ) 

G 

§ 

H 

*( ) 

R = Rr. 



^ ) ( , ) 



Coriolis acceleration 

Acceleration of base with respect to Earth 
Acceleration of base with respect to inertial space 
Indicated position (latitude or longitude) of vehicle 
Indicated position of vehicle at start of problem 
True position of base 
Position correction 

Position correction, at start of problem 

Correction to the indicated vertical 

Correction to indicated vertical at start of problem 

Correction to reference vertical 

Correction to reference vertical at start of problem 

Output signal voltage of device indicated by subscript 

Gravitational field vector 

Gravity field vector 

Angular momentum of gyro unit rotor 

Output current of device indicated by subscript 

Torque on gyro gimbal 

Radius of the Earth, a function of geographic position 
Set value of Earth radius used for Schuler- tuning 

Sensitivity of component or component -group indicated 
by first subscript, with input and output indicated by 
second subscript 
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Indicated vertical 



(vp^ Indicated vertical at start of problem 

V Reference vertical 

r 



True vertical 

(V^Q True vertical at start of problem 

Angular velocity of indicated vertical with respect to 
the Earth 



^(EV)r 

^(EV)t 



W, 

W. 

W, 



IB 

IE 

(IV)i 



cmds 

(U) 

(I) 

'w 

hfi 

Ifi 

Vw 

Vc 

(U)W(g„) 

dir 

fb 

e^ 



Angular velocity of reference vertical with respect to 
the Earth 

Angular velocity of true vertical with respect to the 
Earth 

Angular velocity of base with respect to inertial space 
Angular velocity of Earth with respect to inertial space 

Angular velocity of indicated vertical with respect to 
inertial space 

Controlled member drive system 
Uncertainty in quantity immediately following. 
Inaccuracy in quantity immediately following. 

Water speed compensation current 

High frequency integrator 
Low frequency integrator 
Velocity of craft through the water 

Velocity of water current 
Uncertainty in gyro drift 

Direct channel amplifier 

Feedback amplifier 

Earth rate compensation voltage 



Figure 1-3 represents the functional diagram of the instrumentation of 
a single axis of the complete vertical indicator of Fig. A-1. Figures A -3 
and A -4 are the mathematical functional diagrams of Fig. 1-3 for the un- 
damped and damped systems, respectively. 

The quantity [(C)V]^^ is defined by Equation (A-1); that is, [(C)V]^^ ^ 
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Fig. A -4 Mathematical functional diagram of a basic single -axis damped vertical indicating 
system using the components shown in Fig. 1-3. 



is the angle between the true vertical and the indicated vertical. 
Then, in scalar form and operator notation, 

- W^gyj., 

where 




(A -20) 
(A -21) 



Equation (A -20) corresponds to Equation (A -2) with the coupling term 
(W(Ep)vi X i)) neglected. 

The proposed system consists essentially of an accelerometer, an in- 
tegrator, a single-degree -of -freedom integrating gyro unit, and a drive 
system for the stable platform or controlled member. The accelerometer 
is mounted on, and near the axis of rotation of, the controlled member. It 
will therefore sense the indicated acceleration of the base with respect to 
inertial space, or the projection of the resultant specific inertia reaction 
force on the indicated horizontal plane (in general this force will contain 
both acceleration and gravity effects). 

The development of the system equations is shown in Derivation Sum- 
mary A-I. Equation (A -32) represents the solution of the vertical indicator, 
which in not presented here as an analytical expression but is rather pre- 
sented as a series of figures in Chapter 4, and is further discussed there. 



INSTRUMENTATION OF A SINGLE-AXIS OF THE VERTICAL INDICATOR 
OF FIGURE A-1 (DAMPED). 

Figure A -4, the mathematical block diagram of a damped single-axis 
vertical indicator, is essentially a modification of Fig. A -3, the undamped 
vertical indicator, in which a direct channel and second integration with a 
feedback loop are added in the orientational control-signal -generating loop 
to provide damping and stabilization in the system. The essentials of the 
procedure are as presented in Appendix C. 

The development of the equations for instrumentation of this system 
proceeds in an anaiagous manner to the development in Part H of this ap- 
pendix and is presented briefly in Derivation Summary A -II. 

Equation (A-?8) represents the complete mathematical performance 
equation of the damped vertical indicator. The solution is not presented 
here but is discussed further in Chapter 3. 
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DERIVATION SUMMARY A-I 



The performance equation for the accelerometer unit is (neglecting 
accelerometer dynamics ♦): 

*(au) = S(a.)(a.e)<sC(C)^](t.i) * * ('J)«(av.) • (^^2) 

where e^^ = accelerometer unit output signal , 

S/ w V = accelerometer unit sensitivity for an acceleration 
k TIM » ) input and a voltage output , 

ajg = acceleration of the base with respect to inertial space . 

Since the system intrinsically minimizes [(C)VJ^^ this angle will 
he small enough in general to permit the assumptions : 
sin[(C)T](t.i) = [(C)T](t.i) 

C03 S 1 

In Reference 8 it was shown that, at submarine velocities » 

^ ^^'^(EV)t * (A- 23 ) 

Substituting this in (A^22) gives 
*(au) “ ®(au)(a,e) 

The integrator performance equation is given by 

^(int) ' ®(int)(e,i) p (®(au) * ^'’^®(int) (in)) * ^(int)^; 

* To incorporate Schuler -tuning characteristics, the system is designed 
to have a natural period of approximately 84 minutes. Therefore, the dynamics 
associated with the accelerometer unit and the controlled member drive are 
neglected. 

♦♦ In this scalar equation the centripetal acceleration due to the Earth *s 
daily rotation is included in g rather than The net vertical force is 

effectively g for submarine operation. 
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where : 



^(int) ” current output of integrator « 

Sr. . w \\ ~ seneitivity of integrator for voltage input, current rate 
k n M®, ) outputf 

®(int)(e.i)<"^*(lnt)(ln) = i" l(int) voltage 

input** 

^(int)o ” integrator current at start of problem • 

The primary purpose of the gyro unit is to provide a means of stabi- 
lizing the controlled member against arbitrary motions of the base; however, 
in the present analysis, the gyro unit is simply considered to be part of 
the angular velocity generating system. The input current produces a torque 
on the gyro unit which is balanced by the angular velocity of the controlled 
member with respect to inertial space, thereby keeping the gyro unit on null. 



where 



(tJ)M(gu) = H(U)W^^j ss uncertainty torque on gyro gimbal , 

Sy. w. s= sensitivity of gyro unit torque generator for current 
^ &M * J input and torque output » 

i^ = current compensating for rotation of earth with respect 

to inertial space » 

H = angular momentum of gyro unit » 

** angular velocity of controlled member with respect to 
inertial space . 

W — y 

(IV)i“ (cmds) = angular velocity of indicated vertical with respect 
to inertial space* 



Let 



^(ciDdsHl.W) 



H * 



(A-27) 



* c.f. Chapter 2 - Uncertainties in Components 
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where 



S/ , w. = sensitivity of controlled memher drive 
» ; system for a current input and an angu- 
lar velocity output « 

From the foregoing, 

”(lV)i “ *(S7)i * *IE “ ®(cioas)(i,W)^^(int)* V * 

”lE 

i^ = ^ (I)i^ = compensating current + inaccuracy 

®(cmds)(i,W) in compensating current, (A- '’9) 

*'(EV)i ° ®(cEd8)(l,W)^^(lnt) * 

Substituting the above relationships into (A-30) , we obtain 
^(EV)^ " ^ ^ ''(SV)t " “ ^(EV)i^ 

= ®(cMd8)(l,W) [®(au)(a.e)(eC(°)V](t.i) * V 

* (^)«(au) * W)«(lnt)(in)] * ^c^s)(l.W) ^ 

^(cmds)(i,w) ^(int)o * (A-.3I) 

Equation (A-31) corresponds to Eq. (A-4) . Substituting E q, (A^31) into 
Eq. (A-20), differentiating the result and setting 

®(au)(a,e) ®(int)(e.i) ®(ouds)(l.¥) = T we 

obtain the differential performance equation 






®e®(au)(a,e) 






-P[^c»de)(i.W) W)'V) 



(>-33) 
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Derivation Summary A- I I 



The performance equation of the accelerometer unit is unchanged "by 
the modifications cmd is: 

®(au) ° V(ST)^] 

However, the output of the integrator is now: 



" ®(ni) * ^”^®(hfi)(in)) * ^(int)o ’ 



where 



®(dlr) ° ®(dir)(e,e) ®(au) * ^°^®(dir) ’ 



®(int) * ^ (®(au) ■*■ ^”^®(lfi)(in))* ®(int)o ’ 

*(cp) ” voltage compensating for ship*8 speed through the water, 

®(fh) " ^(fb)(i.e) ^int ^^^®(fb) » 



” performance function of hi^ frequency integrator, 
” performance function of low frequency integrator . 



The torque “balance eq\uation for the gyro \uiit remains unchanged. 
* S(tg)(i.M) (^int ♦ V ' H«(cada) ’ ' 

, V, obtain 

’'(IV)i “ ’'(E7)i * *'(11!) “ ''(lomde) ’ 



(A-26) 

(A-27) 

(A-28) 
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"" S 



"IE 
(cmde) (i,w) 



+ (I)i . 

' ' w » 



• (A- 29 ) 



«(OT)i = S(c„de)(i.W) <^(int) ♦ ^ <^)«(gu) ’ 

Sutatituting the above relationships into Eq.. (A-30) gives 
* ®(hfi)(e,i) ®(fl>)(i,e)^ ("(EV)^ - P 
° ®(hfi)(e,i) ®(cind8)(i,W)( ®(dir) * ®(int) * ®(cp) ^^^®(cp) 

- («)®(ft) * (“)®(hfl)(in)V((«V) * ^cbcIbHIW)) 

x(P * ®(hfi)(e,i) ®(ft.)(i,e)) 

“ ®(hfi)(e,i) ®(cmde)(i,W)[( * ®^"(ET)^) 

q 4. ^(lfi)(e,e) ""(au)(a.e)\ 

V (dir)(e,e) ^(au)(a,e) ^ p / 

^ ^ ^(lfi)(e.e) \ f... . ^(lfi)(e.e) 

V (dir) (e,e) p ®(au) p '®(lfi)(in) 

^^^®(dir) ■*■ ®(cp) ^^^®(cp) ■ ^^^®(hfi)(in)*" ®(int)o] 

■"(P * ®(hfl)(e,i) ®(ft>) (i,e)X“^''(gu) ■" ®(eM9)(i,W)<^^^)- 

(A- 3 i() 

Setting ^(cmde) (i,W) ^(au)(a,e) ^(dir)(e,e) “ R 

® (A^ 35 ) 

Schuler tuning and differentiating (A-34) and after re-arranging, we get 

(P® *P* ^hfiXe.i) ®(n.)(i,e) *P ^ 

r-A H /'c • <; ” ^(Ifi).( e.elo 

= [p-(l - r p (S(hft)(e.i) ®(n,)(i.e) - 
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(U)e 

^8 ^(au)(a,e) 



i& u) _ / ^(lfl)(e .e)\ _ . . , 

(aulfa.e)''^ ^ (hfi)(e.i) (ci»ie)(l.w)( <”^*(dir) 



’(dlr)(e,e)' 

* «(cp) * (I)«(cp) - * (^)»(hfi)(in)) 

• ®(hfi)(e,i) ®(cDda)(i.w) ®(lfi)(e,e)^’’^®(lfi)(ln) 

- (p® ♦ P S(hfi)(e,i) ®(fi,)(i,e)) - ®(ci.d»)(i,w) 



With = E very nearly, Wriglsy(^) has shown that the tem^l - ^'^will be 

negligitly email under conditions which will he encountered in service. In 
Appendix (C), it is shown that it is possible to minimize the term in 

P[s(hfi)(e,i) 6(fb)(l,e) - "(M), 6anP^”« 

ratio, 

«(E7), = r*r- 

where = component of craft velocity with respect to the water perpen- 
dicular to input axis of system., 

and 7^ = component of current velocity with respect to the earth perpen- 

dicular to input axis of system. 



/s • S - 

/®(hfi)(e,i)®(fb)(i,e) S, w Av 

Now setting e, , = s r— s LdirJie^ j 

cep; \ ®(hfi)(e,l) ®(cmdB)(i,w) /® 



(A-37) 



the terms in (A-36) leaving only a term 

P(6(hfi)(e,i) 6(fb)(i,e) - r 
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will In general "be negligil)le, 'because '7^, the rate of change of current 
velocity, is in itself generally small, and, in addition, because the coef- 
ficient of is optimized (Chapter 4). Note that the effect of any errors 
in the measurement of (V^) are considered In the above system equations in 
the term (I)e^^p^. 

Equation (A^36) represents the differential performance eq^jation for the 
system. .If now the Laplace transform of the performance equation is taken, 
the initial conditions appear as follows where now p^Laplace operator. 



S s 



(p® + P® ®(fT>)(i,e) * P i" * S 



'a ®(dlr)(e,e) 






(C)7] 



(t.i) 



- ^ l®(hfi)(e.i) ®(fb)(l.e) - js 

^(dir)(e,e) ^ ^s ^(au)(a,e) 

■I* ®(hfl)(e,i) ®(cind9)(i,w) [^°^®(dir) 



* ^°^®(hfl)(ln)j "®(hfi)(e,i) ®(cinde) (l,v») ® (Ifl) (e,e) ^®^® (Ifl) (In) 



”(P ^®(hfi) (e,i) ®(fb)(i,e)) ® (cmds) (l,w) 

X (p® + P 

®(hfi)(e,i) ®(ft)(i,.)) + (P® P ®(hfl)(.,i) ®(fb)(i,.) 



8 

(U)e 



(au)o 



R S 



(au)(a,e) 



+ S 



(hfi)(e,i) ^(cmds) (i,w) ^ ^^^®(dir)o ^^^®(cp)o 
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* ®(hfl)(e,i) ®(fb)(i,.)) ^”^*'(su)o * ®(cmae)(i,w) 
* ®(cisdB)(i,w) * ®(hfl)(a,i) ®(fli)(i,e) ) • 



(A-38) 
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APPENDIX B 



POSITION INDICATION BY OPEN-CHAIN INTEGRATION 
OF THE INDICATED VERTICAL 



Figures B-1 and B-2 are the mathematical functional diagrams of 
the proposed position indicators corresponding to the damped and un- 
damped vertical indicators respectively, as described in Appendix A. 

Derivation Summary B-I presents the mathematical derivation of 
position indication as derived for both systems. This derivation has 
the same form whether the vertical indicator is damped or undamped. 
This fact is readily apparent when it is recalled that the relationship 
between the angular velocity of the controlled member and the inte- 
grator current is the same regardless of damping. The relationship 
is given by Equation (A -30) . 

Equation (B-6) shows that the correction to the indicated position 
is the sum of the correction to the indicated vertical, initial conditions 
in position and correction to the vertical, together with the time integral 
of certain uncertainties and inaccuracies. 

Recalling that the transfer function has been defined as the response 

of the system to the forcing function we may write |f(w, t)| for 

the open-chain integrator as 1 - 1 - (1 - The mean-square er- 

I I 

ror may then be evaluated quite simply after the method applied to the 
undamped closed-chain vertical indicator in Chapter 3. That is, the 
plotof|^(l-eW)i2 may be constructed and the integral of Equation 
(2-11) apjAied. 



^(output) 



t = t, 




G(i„)(w) dw, 



( 2 - 11 ) 



where represents the power spectral density of each of the uncer- 
tainties shown in Equation (B-5) taken in turn. 
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(int) 
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Fig. B-1 Mathematical functional diagram of an undamped system for indicating position by 
means of an open-chain integration of the angular velocity of the indicated vertical. 



(do) 
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Fig. B-2 Mathematical functional diagram of a damped system for indicating position by means 
of an open-chain integration of the angular velocity of the indicated vertical. 



Derivation Summary B-I 



Position Indication "by Open- Chain Single Integration 
of the Indicated Vertical 



In both the systems of Figs. B-1 and B-2, position information is ob- 
tained from a single integration of the current input to the gyro unit. 

Writing the differential performance equation for the open chain integration: 

where ~ uncertainty in output of the position indicating 

' '' ' integrator reduced to an equivalent signal input., 

(P^)^ = initial position. 

From^q. (A-30), and noting that (A-20), 

- 'C<'>’)(..D - <"»(„))- • 

(B-2) 

Substituting Bq. (B-2) into Bi. (B-1) and arbitrarily setting 
S(cMs)(i,w) =®(Pi)(iP)- 

^1 = ? («(EV)t - P C<‘=5’J(t,i) - * <^^^(Pi)(in) - 

-S(Pl)(i,P) (I)V)*(Pi), . (B-3) 

noting that i = P^ -(P^)^ , we write 



i(p[(0)V](,^l)) =[(0)V](^^i) -C(0)V](t^i)o 
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(B-») 



"?«°)‘(Pi)(in)- (°)^gu)-^Pi)(i.P) • 

Further tdefining [(C)Pj^^ i) “ * 

C^°)P^(t,i)o “ "^1^0 • 

C(c)p](,^i) = [(c)P](t,i), * [(C)V](,^,) - C(c)7](,,i), 

■ I ''(gu) - ®(Pi)(l,P)(^>l^) • 

Writing this in the more conventional form we get 

- ^o((")^(Pl)(in) - ("^■-'(gu) - S(Pi)(iP) 



(B-5) 



(B-6) 
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APPENDIX C 



OPTIMIZATION OF SYSTEM PARAMETERS 



In any dynamical system, the engineer is forced to a compromise be- 
tween a system with minimum forced or steady-state errors and a satis- 
factory transient response. 

Given a performance fimction of the form 



[PF] 



CqP +C^p + — 
^^in’^out^ (1 + Sq) p*' + Sj^p”” ^ + - 




(P) , (C-1) 



Blasingame^^) has concluded that an optimum practical system may be 
achieved with a performance function whose characteristic equation is 
a cubic; and furthermore, such a system may be achieved by a network 
involving only integrators and summing devices, as shown in Fig. C-1, 
where the C’s and S’s may have all values, positive or negative, including 
zero. 

Such an optimum system has been applied to the undamped vertical 
indicator of Fig. A- 3 to obtain the damped vertical indicating system of 
Fig. A-4. 
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^(out) 



(t) 



-©-i 



. coefficient multipliers 

Fig. C-1 Generalized mechanization computer. 

The performance equation for the vertical indicator of Fig. A-4 for 
an input in the angular velocity of the true vertical may be written 



r , D^ri-S 

[(C)V]jy^ _ P ^(Vi)(a,W)^>-^PVj Tn 



W 



(EV). 



p2 + 2(DR)W„3P + w2^ p+^ 



, (C-2) 



P ®(Vi)(a,W)^ + P^ T„ J . (t, J) 



P^+-^P^ + rP +■“ R 

Td R R 



where 



^(Vi)(a,W) = ^^'^(a,e)^(dir)(e,e)^(hfi)(e,i)^(cmds)(i,W) 



2(DR)W, 



ns+'TO: = S 



^(lfl)(e,e) 



(dir)(e,e) 



(C-4) 
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1 


'^ns 1 


'^n 


w2 (CT), 
'^nE * 


(ct)d = 


2(DR)W„^ 





®(hfi)(e,i)®(fb)(i,e) 



(C-5) 



(C-6) 



" nE ■ " ns 

Wrs = natural angular frequency of the system, 

WnE = earth’s natural angular frequency, 

(CT)j^ = system characteristic time. 

2 

The p term in the numerator of Equation (C-2) has the significance 
of an acceleration in the angular velocity of the true vertical and gives 
rise to errors known as jerk- errors in the system. The p term in the 
numerator of Equation (C-2) has the significance of a rate of change in 
the angular velocity of the true vertical and gives rise to errors in the 
system known as velocity errors. 

The physical significance of the two changes associated with the angu- 
lar velocity of the true vertical may be seen if it is recalled that movement 
of a craft over the surface of the earth produces a change in the direction 
of the true vertical. This change in direction or angular velocity of the 
true vertical may be represented as a vector at right angles to 

the radius vector of the earth and the velocity vector of the craft, and of 
I V. _ 

^ , where is the velocity of the craft and is the 
E 

radius of the earth. 

Thus, a constant linear acceleration of the craft over the earth re- 
sults in a steady rate of change in the magnitude of the vector , 

and is physically associated with the p term in the numerator of Equa- 
tion (C-2). However, maneuvering accelerations normally encoimtered 



magnitude 
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on shipboard, say in a change of speed from 5 to 15 knots, involve a rapid 
change of acceleration from zero to a maximum of about 5 knots per min- 
ute, with another rapid drop in acceleration as the new speed is approached. 

2 

These changes in the rate of acceleration are associated with the p term 
in the numerator of Equation (C-2), and numerically are considerably 
greater than the maximum accelerations involved. Similarly, a change of 
course involves a change in the magnitude of the components of 
along each of the two horizontal axes of the vertical indicator that is sta- 
bilized in azimuth. Thus, a course change also gives rise to both jerk 
and acceleration errors, in which the rate of change of acceleration will 
again exceed the acceleration for a normal tactical maneuver. 

An examination of Equation (C-2) shows that it should be possible to 
eliminate one of the terms in the numerator by a proper choice of sensi- 
tivities, and then to minimize the effects of the remaining term. A study 
of acceleration and jerk effects likely to be encountered on shipboard in- 
dicates that it may be desirable to proceed on the basis of eliminating the 
jerk- error term and minimizing the acceleration-error term, as has been 
done in the following sections. 

With the jerk-error term eliminated by setting W) " R ^ 

Equation (C-2) may be written 



^(EV)^ 





P^2(DR)W„^ + w2^^(P+^) 



(C-7) 



Letting 

pjWf = j(FR)W^g , 

= angular frequency of forcing function , 
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(FR) s 



[<C)V] 



M = 



Wf 




^ ^nE ’ 




^ns 




' ^nE ’ 




1 (C-7), we 


obtain 


i(FR)W„j, 


4(DR)2(FR)2^ + (i - (FR)2„ 


(CT)j3 


1 - {FR)\ 



^(EV)j W^„^(FR)2 



f . im . i . j ? (D R)irR). / , /yp^w frx^ 4.1 
(CT)p (FR)„ + ^ 

(C-8) 

Hydlnger and Webster have shown that the forcing frequencies in 
W(Ev) likely to be encountered on shipboard will be essentially of zero 
frequency or such that (FR) >10. 

Optimizing on the basis of (FR) > 10, it may be shown that the following 
approximations hold: 



a. (FR)W„g(CT)j 3 >> 1 



b. 



1 . iER)lU> 2fDRV-iEB 
(FR)2r ^ W 



as) 

n 



c. JEEf » 1 



(FR) 



By making use of these approximations, Equation (C-8) is reduced to 
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[(c)v]ft.n 



4(DR)^(FR)^„+ (l-(FR)^„) 



2n 



1 - (FR)** 



(FR)W„e 



w2„e(FR)2„ 



. (g R)-- 

(FR)2„ 



j (FR)W„e(CT)j3 



(C-9) 



<CT)dW"„e 



(FR)-" 



4( DR)^(FR)^„+ (l - (FR)^J 

i-(frA 



(C-10) 



Substituting 



(CT) - _ 2(DR)(FR)„ 

' <E-<s ’ W„eU-(FR)2] 

and rearranging, we obtain 

[(C) V] 



(C-6) 



W(EV)f ^ 2 (FR)V„e (DR)(FR)„ - ^ ^"(DR) 

(c-11) 

Assuming a fixed (DR) and taking the partial derivative of Equation 
(C-11) with respect to (FR)„, we have the necessary condition for a 
minimum in I to exist at a given damping ratio. 



A(C)V]ftiA ^ 

1 to exist at a given damping ratic 

\^(EV)^ / 



I - 0 



2(FR)'»„j 



4( FR)^(DR)^-1-2(FR)^ + 3(FR)^ 
(FR)^(DR) 



(C-12) 
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Solving Equation (C-12) for (DR) as a function of (FR)^^, we have 



(DR)‘‘ = 



2 1 + 2(FR)2 -3(FR)t 



4(FR)‘ 



(C-13) 



Figure C-2 is a plot of (DR) vs. (FR)^^ from Equation (C-13). With 
(FR)j^ = 1, (DR) = 0, and with these conditions in Equation (C-13), we 



find 



[(C)V] 



W 



lUi 



(EV). 



= 0. The conditions are essentially those of an undamped 



Schuler- tunei inertial system, which Wrigley^^^^ has shown to be insensi- 
tive to a sinusoidal forcing function with an 84.4-minute period correspond- 
ing to (FR)^^ = 1. 

However, a condition of zero damping and (FR)^^ = 1 is neither 
attainable nor desirable in practice for an inertial system operating over 
an extended period. Hence, choosing (DR) =0.1, 0.3, and 0.5 for a com- 
parative study, and utilizing the relationships of Equations (C-16) through 
(C-18) as shown below, we arrive at Table C-1, which is a tabulation of 
the system parameters for the vertical indicator of Fig. A- 4, as obtained 
when optimum values of (FR)j^ are chosen for the various damping ratios. 
The resulting characteristic equations are also shown below. 

The characteristic equation of (C-2), which is in the form of a cubic, 
can also be expressed as the product of a second-order term and a first- 
order term as follows: 



(P^ + 32?^ + p + = (p^ + Ap + B)(p + D) ; 

(C-14) 

also as 

(p2 + AP + B)(p+D) = (P^ + 2(DR)W„gP + w2g)(p + -^^) . 

(C-15) 
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o 




p: 

A 
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(FR)„ 

Fig. C-2 Optimum (DR) as a function of (FR) 



Equating coefficients of like terms, we get 



ag = (A + D) = 2(DR)Wj^g + 



1 



ns + (CT)jj * 



(C-16) 



= (B + AD) = + 



2(DR)W„, 
(CT)p > 



(C-17) 



ap = (BD) = , 




(C-18) 



(CT)jj 



2(DR)W„^ 2(DR)(FR)„ 



(C-6) 




Figures C-3 and C-4 are the resulting steady -state responses of 
the optimized vertical indicator to sinusoidal forcing motions of the 
true vertical. It will be noted that the response to accelerations is 
three orders of magnitude greater than the response to angular velo- 
cities at (FR) = 1; however, the response to accelerations is decreas- 
ing at 18 decibels per octave, whereas the response to velocities is de- 
creasing at 12 decibels per octave. Consequently, at (FR) - 100, the 
responses to both inputs are of equal magnitude, while at higher fre- 
quencies the response to accelerations becomes increasingly lower 
than the response to velocity inputs. 
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JM/suo|poJ!iijUJ 



SUDlpOJJIlim ^’‘*^[A0)] 
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Fig. C-4 Steady-state response of vertical indicator to sinusoidal angular acceleration 

of true vertical. 



APPENDIX D 



METHOD USED IN THE DETERMINATION OF SYSTEM 
PERFORMANCE CHARACTERISTICS 



Consider a generalized performance equation of the form 



(p3 + E2p2 + a^p + a^,) q(p)„^t = <l(P) ^ 

+p(p + a^)q(p)(i„^) - - 1(in3)(°+) ■ ^ 

■ 

where 

p = Laplace transform variable , 
q(p)|n ^ Laplace transform of q(t)jj^ , 
q(in)(o+) = initial value of q^^^j at t = (o+), and 
a^, a^, a2> . . are numerical constants. 



The complete response of the differential equation may be found as the 
sum of the responses to the various terms on the right-hand side of the 
equation. Thus, the generalized performance equation may be broken down 
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into a series of performance equations with but a single term on the right, 
such that q(p)Q^t may be represented as follows: 

Setting 

[P] = ( p^ + a2P^ + a^^p + a^) , the characteristic equation, 

we obtain 

1<P)(out^) = > (D-2) 

<l(P)(out2) = 

'J^^outg) = [pf 

q(p)(out^) and so forth. 



Expressions (D-2), (D-3), and (D-4) show that the generalized system 
may be considered as made up of a number of independent sub-systems 
with performance equations given by 

1 p . 

3 2 >32 > 3 2 > 

p +a2P +ajp + a^ p +a2P +a^p + a^ p +a2P +a^p+ a^^ 



respectively. Thus, the weighting function for the generalized system 
for an input in q(t)^^ j may be determined quite simply by setting the sub- 
system of q(tL X into an analog computer and taking the time response to 
' n' 

a imit impulse. 

It should be noted here that the terms which have to do with the 
initial conditions are in the form of a product of a performance equation 
and a constant, and that the inverse Laplace transform of a constant has 
the form of an impulse; hence the time response to initial conditions has 



87 



the form of a weighting function. Similarly, if some probable 

form that may be predicted on the basis of its physical nature, the probable 
response of the physical system to may be found. In particular, 

if the system inputs and outputs are represented as in the mathematical 
block diagram of (D-5), 



then 



where 



aje 



jWjt 



(PF) system 



ajF(Wf,t)e' 



jWft 



(D-5) 



q(t) 



(out) 



= ajF(Wf,t)e’ 



jWft 



af|F(Wj,t) e 



J(Wft + A^^) 



dyn' 



(PF)gyg = system performance equation of form of (D-2) etc.. 



F(Wf,t) 

■^dyn 



transient transfer function for frequency w- (c .f . 
Chapter 2), 

d 3 mamic phase angle associated with F(w,t) ; 



and, if it is desired to find the absolute value of I(w£,t) as a function of 
time for a given frequency w^, then the output of the system to the forcing 
function 

jw.t 

aj e = a^(cos W£t 4 jsin Wjt ) , (D-6) 



for which the system response may be written 

= aj|F(Wj,t)|cos (Wjt+ A^yjj) + jaj|F(Wf,t)| sin (Wjt + A^yjj) , 

(D-7) 

may be represented on the complex plane as in Fig. D-1. 
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Fig. D-1 

Complex plane representation of F(wpt)e^^f^. 
From Fig. (D-1), it is apparent that 



I f f I I \/cos2(Wft H- A^y„) + sin2(Wft + , 



Solving for | F(wj,t|, we get 

= I ^^^^out 

ti aj 






t = t, 



(D-8) 



(D-9) 



The mathematical relationships of the foregoing paragraphs suggest 
that |F(w£,t)| may be obtained quite simply by use of a high speed analog 
computer utilizing timing marks to indicate time and clamping to ensure 
zero initial conditions at the start of each computing cycle. 

If the input to the system is chosen as a sinusoidal input of unit 
amplitude such as cos w^t, the output will be obtained as a function of 
time such that 




(D-10) 



If then the output (D-10) is put through an integrator with the performance 
function w^p , the result will be 
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(D-11) 




COS Wjt dt 

0 

= |F(Wf,t)| sin (Wft + AjjyJ. 

If the output (D-10) is now taken to the x^axis input of an oscilloscope 
and the output (D-11) is taken to the y-axis input of the oscilloscope, 
the picture on the oscilloscope will be the time varying presentation 
of Fig. (D-1), where the radius from the center of the diagram to the 
resulting curve is representative of the sum in Equation (D-7), which 
has been shown in Equation (D -9) to be equivalent to 



or 



I F(Wf,t) I 



q(t)out 



I F(Wj,t) I 



‘iWoutl 



= radius of oscilloscope curve at 



time (t), for a sinusoidal input amplitude. 

The function | F(w,t)|, which was derived in Chapter 2, may thus be 
determined by applying unit forcing functions of various frequencies 
over a frequency band for which the response is desired and taking the 
amplitude of the radius of the oscilloscope presentation as a function of 
time. 

It is apparent from Fig. D-1 that knowing the frequency and the time 
for any point on the curve, and having a base line established by some 
means, the dynamic phase angle may be easily determined from the 
relationship 



A(dyn) =A - V' 
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